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Abstract 

We report a 13 C-NMR study on the organic conductor <9-(BEDT-TTF) 2 l3 [BEDT-TTF: 
bisethylenedithio-tetrathiafulvalene] at ambient pressure. The lattice symmetry is consistent with 
the pseudo-orthorhombic structure by the orientation dependence of the NMR spectra, in which 
all molecules in the unit cell are crystallographically equivalent. There is no signature of charge 
disproportionation as observed in most of the #-type BEDT-TTF salts. The analysis of NMR 
Knight shift, K, and the nuclear spin-lattice relaxation rate, 1/Ti, revealed the degree of electron 
correlation, evaluated by the Korringa ratio (oc 1/T\TK 2 ), is in an intermediate regime. However, 
NMR relaxation rate is enhanced above 200 K, which is a possible indication that the system enters 
into a quantum critical regime of charge-order fluctuations as suggested theoretically. 

PACS numbers: 
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FIG. 1: (Color online) (a) Side view of the pseudo-orthorhombic crystal structure of 6-I3. (b) 
Top view of the schematic arrangement of BEDT-TTF molecules in the conducting ac plane, (c) 
Molecular principal axes of a BEDT-TTF molecule. Arrows indicate the positions of 13 C atoms. 
In figures (a) and (b), two molecules which become unequal in their geometric arrangement with 
respect to the field direction of H / / be and H / / ac are represented by dotted and dashed ellipses 
in (a) and (b), respectively. 

I. INTRODUCTION 

The layered organic compounds (BEDT-TTF) 2 X exhibit fascinating electronic phases 
with superconductivity-^, Mott localization^, spin liquid^, charge ordering^, massless Dirac 
fermions^, and so on, where X is monovalent anion and BEDT-TTF is bisethylenedithio- 
tetrathiafulvalene. The system is constructed fom alternate stacking of 2D conducting layers 
of the donor molecule BEDT-TTF and insulating layers of anion X [Fig. Ufa)]. The anion 
extracts an electron from two BEDT-TTF molecules, i.e., (BEDT-TTF) +05 : X~, which 
results in a quarter-filled hole band through inter-molecular overlaps of the molecular orbital. 

The title compound, 9- (BEDT-TTF) 2I3 (abbreviated to 6-I3), is known as a typical 2D 
metal with highly symmetric zigzag alignment of BEDT-TTF molecules in the conducting 
layer [Fig. H^b)] 8 . It is metallic in the whole temperature range 8 , and the Fermi surface has 
been confirmed by quantum oscillation^— and optical^^experiments. Magnetic suscepti- 
bility exhibits weak temperature dependences^, characteristic of the Pauli paramagnetism. 
Moreover, superconductivity is observed at Tc of ca. 3.6 K-^&. 

Meanwhile, the isostructural materials #-(BEDT-TTF) 2 MZn(SCN) 4 [M = Rb and Cs] 
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(abbreviated to 0-MZn), undergo metal-to-insulator (MI) transitions at Tui = 195K (6- 
RbZn) and 20K ((9-CsZn) 4 accompanied by a spontaneous charge ordering 5 (CO) and a 
glassy charge disproportionate!^, respectively. Even in the conducting state above T MI , 
where resistivity is temperature insensitive 4 -, 13 C-NMR line shows a gradual broadening 5 -^ 
toward Tmi, which is a manifestation of charge fluctuations. Thus, the electron correlation 
plays a decisive role in the electronic phases in #-MZn. Despite the close similarity in the 
molecular arrangements between (9-MZn and 6-I3, their contrasting ground states demon- 
strate a variation in the effects of electron correlations depending on the bandwidth, as 
pointed out by Mori et al.—. 

Application of hydrostatic pressures drastically alters the electronic state in (9-I 3 . Pressure 
dependence of the room-temperature resistivity shows a discontinuous jump at around 0.5 
GPa 20 . Under pressures above this critical value, 8-I3 shows weak temperature dependence in 
resistivity- 1 , a large increase in Hall coefficient^, and Independence in nuclear spin-lattice 
relaxation rate, 1/T 1; below 20 K— with decreasing temperature, indicating the presence of 
a linearly dispersive energy band near the Fermi level^ 1 ^ 4 -. Notice that the same properties 
are also observed in a-(BEDT-TTF)2i3 (a-13) at high pressures approximately above 1.5 
GPa- 1 ^ 5 -, where a tilted massless Dirac cone is predicted by band-structure calculations^ 1 ^. 

includes three nonequivalent molecules in the unit cell with the zigzag molecular ar- 
rangement as in the 9 type^. It is suggested that this a-type local site symmetry plays 
a significant role in the realization of the tilted Dirac cone in a-I^-^^. Similar argument 
is expected in #-I 3 under pressure, though the crystal structure at above 0.5 GPa is not 
yet determined. Moreover, x-ray diffraction experiments at ambient pressure suggest two 
possible types of crystal structures in (9-D/, i.e., pseudo-orthorhombic lattice (space group 
-Pnma) and monoclinic lattice (P2i/c). This open issue should be explored from different 
experimental probes. In order to reveal the electronic states in 8-I3 comprehensively, it is 
hence important to make a microscopic characterization of the ambient-pressure metallic 
state and the local site symmetry. 

As seen above, (9-13 is a unique material to connect the physics of correlated electrons 
and Dirac electrons, which are both among the intensively studied issues in condensed 
matter physics. In the present work, we aim at elucidating the electronic properties in 
6-I3 at ambient pressure by means of 13 C-NMR experiments for the first time. This is 
expected to provide a basis for understanding the conducting phase in the proximity of 
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CO phase and the massless Dirac fermions in organic conductors. First, we measured the 
orientation dependence of the 13 C-NMR spectra against applied field to evaluate the local- 
site symmetry and the hyperfine-shift tensor at the central 13 C sites in BEDT-TTF molecules 
[Fig.[TJc)]. Because the hyperfine interaction between electron and 13 C-nuclear spins is highly 
anisotropic in (BEDT-TTF) 2 X compounds^, the hyperfine-shift tensor 5 = (S xx , S yy , 5 ZZ ) at 
the 13 C sites is thereby determined, where x, y, and z represent the principal axes of the 
molecule [see Fig. QJc)]. Then, we measured the Knight shift, K, and the nuclear spin-lattice 
relaxation rate, 1/Ti, and determined local spin susceptibility and electron correlations in a 
quantitative manner, using the hyperfine-shift tensor determined. Based on the experimental 
data and analysis, we discuss the nature of electronic state in 9-l 3 at ambient pressure. 

II. EXPERIMENTAL 

Single crystal of 8-I3 was prepared by the conventional electrochemical method. For the 
13 C-NMR measurements, the central carbon sites in BEDT-TTF molecules were selectively 
enriched by 13 C atoms (nuclear spin / = 1/2, 7 n = 10.7054 MHz/T) with 99 % concen- 
tration [Fig. Hfc)]. All NMR measurements were performed for a single crystal at ambient 
pressure in a magnetic field H of approximately 6.00 T, which is rotated within the con- 
ducting ac plane and the be plane [see Fig. [TJ. NMR signals were obtained through the 
fast Fourier transformation of the so-called solid-echo signals, and the 13 C-resonance fre- 
quency of tetramethylsilane [(CH 3 ) 4 Si, TMS] was used as the origin of the NMR shift. The 
nuclear spin-lattice relaxation rate, 1/Ti, was determined by the standard saturation and 
recovery method. The relaxation curves of nuclear magnetization were well fitted to single 
exponential functions over a decade. 

III. RESULTS AND DISCUSSIONS 

A. NMR spectra, local-site symmetry, and hyperfine-shift tensors 

The local-site symmetry at the central 13 C positions of BEDT-TTF molecules [Fig. QJc)] 
is examined by the angular profile of 13 C-NMR spectra, which were measured with changing 
(i) the angle 9 between the crystal a axis and the applied field H in the ac conducting 
plane [Fig. d^b)], and (ii) the angle ip between the crystal b axis and H in the be plane 
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[Fig. Ha)]. Figures [2^a) and|2]^b) show the typical angular dependence of 13 C-NMR spectra 
under H applied within ac plane (at 40 K) and under H applied within be plane (at 100 
K), respectively. Under all field orientations, the NMR spectra consist approximately of 
four lines at all measured temperatures. Two pairs of lines in Fig. |2]^a) exhibit out-of-phase 
angular dependence against 8, whereas an in-phase variation is seen against ip in Fig. [2]^b). 

As we mentioned in Section I, x-ray diffraction experiments suggest two types of crystal 
structures in #-I 3 ; one is the pseudo-orthorhombic lattice (with space group P nma ) and the 
other is the monoclinic one (with space group P2i/c) [8]. In the pseudo-orthorhombic 
symmetry, the unit cell contains two conducting layers and all molecules in the unit cell are 
crystallographically equivalent. Moreover, the inversion center locates exactly in-between 
the two adjacent 13 C atoms in a BEDT-TTF. In the monoclinic lattice, three molecules in 
each layer of the double-decker unit cell are crystallographically nonequivalent and only two 
of them have the inversion symmetry in the molecular center as in a-l^. In the 13 C-NMR 
spectra, each nonequivalent molecule gives distinct resonance lines. Furthermore, due to 
the nuclear dipole interaction between the adjacent 13 C nuclei, each resonance line from 
a molecule splits into a doublet or a quartet depending on whether the molecule has the 
inversion symmetry or not^I. In case of the pseudo-orthorhombic lattice, when H is applied 
within ac plane, the two layers in the unit cell are equivalent while the two differently 
oriented molecules in a layer become unequal with respect to the field H. On the contrary, 
when H is applied within be plane, the two adjacent layers are unequal whereas the two 
differently oriented molecules in a layer become equivalent with respect to H. Thus, only 
two molecules with inversion centers are unequal against magnetic field H in both field 
geometries and four lines (two pairs of doublets) are expected for the pseudo-orthorhombic 
lattice. In case of the monoclinic lattice, however, three molecules (one without and two 
with inversion centers) are unequal in a layer against an arbitrary field direction. Thus, eight 
lines (one quartet and two doublets) are expected under the ac-plan-field geometry and they 
should be doubled under the 6c-plane-field one, where the adjacent layers are unequal. It 
is obvious that the observed four-line structures in Figs. [2]^a) and Mjo) are fully consistent 
with the pseudo-orthorhombic lattice. A small peak splitting in several spectra, e.g., at 9 
= 20°, 30°, 40°, 50°, and 60° in Fig. [2](a), is attributable to a slight misalignment of the 
rotating plane of H from ac plane. Even if the split lines are counted as independent lines, 
the number of lines is much less than what is expected in the monoclinic lattice. 
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FIG. 2: (Color online) Orientation dependence of the 13 C-NMR spectra (a) under H within ac 
plane (40 K) and (b) under H within be plane (100 K). The 9 and ip are field angles from a and 
b axes, respectively. Symbols represent different Pake doublets stemming from the two differently 
oriented molecules in the double-decker unit cell (See Fig. 1). 
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This becomes clearer when we look into the angular dependence of the NMR line structure 
in more detail. Based on the pseudo-orthorhombic symmetry, we assigned the observed four 
lines into two Pake doublets, as denoted by circles and triangles in Figs. ^ a) and^b), 
stemming from the two differently oriented molecules in the unit cell [dotted and dashed 
ellipses in Fig. [TJb) for the lines in Fig. I2^a), and those in Fig. [TJa) for the lines in Fig. 

BD>)1- 

Figures Eta) and [3](b) show the splitting width, d, of each doublet. The solid and 
dashed curves represent the calculated angular dependence of d, which is expressed as 
d = (3/2r 3 )7^(l — 3cos 2 #) 29 , where 7 n = 10.7054MHz/T is the gyromagnetic ratio of 
13 C nucleus, 2irh is the Planck constant, r denotes the distance between the adjacent 13 C 
nuclei, and 9 is the angle between H and the molecular x axis [see Fig. H^c)]. We used the 
crystal structure reported by Kobayashi et al.— and optimal r = 0.135 nrn^i. The calculated 
curves are in good agreement with the observed splitting width. 

Figures Eta) andBJb) depict the angular dependence of the 13 C-NMR shift given by the 
midpoint of each Pake doublet under H//ac and H/ /be at 100 K, respectively (the same 
symbols as in Figs. [2] and [3] are used). We performed a least square fit of sinusoidal curves to 
these data sets in Figs. Eta) and Etb) simultaneously with the three principal values of the 
hyperfine-shift tensor as fitting parameters based on the crystal structure 3 - 3 -. The resulting 
fitting curves are shown by solid and dashed curves in Figs. Et a ) andHJjb), where the principal 
values of the hyperfine-shift tensor are uniquely determined as (S xx , 6 yy , 5 ZZ ) = (67, 134, 
300) (in ppm). Again, the agreement between the calculated curves and the experimental 
data are fairly well. From all these arguments, it is concluded that all of the BEDT-TTF 
molecules in 8-l 3 are electronically equivalent, being consistent with the local symmetry 
of the pseudo-orthorhombic lattice. This agrees well with the quantum oscillation 9-14 and 
optical reflectance^ 1 ^ experiments both of which confirmed a round and large Fermi surface 
predicted by the band-structure calculation with the pseudo-orthorhombic symmetry^. It 
is noted that, if the three crystallographically non-equivalent molecules in the monoclinic 
structure are not electronically distinctive, this case is not ruled out. 

It is also noteworthy to mention that each line in the NMR spectra is much narrower 
(~ 1kHz) in the whole temperature range than the cases in #-RbZn and #-CsZn which 
exhibit NMR line splitting or broadening over several kHz or more in association with 
the charge ordering and its glassy freezing, respectively*^. The present spectral feature 
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FIG. 3: (Color online) Angular dependence of the nuclear dipole splitting width d at the two 
differently oriented molecules in the unit cell (denoted by symbols) under (a) H//ac (40K) and (b) 
H/ /be (100K), extracted from Figs. 2(a) and 2(b), respectively. The same symbols are used as in 
Fig. 2. Calculated angular dependence is shown by solid and dashed curves based on the crystal 
structure reported by Kobayashi et aZ. 33 and the optimal 13 C= 13 C bond length of r = 0.135 nm^i. 

clearly indicates that 8-I3 at ambient pressure is free from these instabilities of the charge 
organization. This is consistent with the sharp line width observed in Raman scattering 
experiments^. 

As mentioned in Section I, the metallic phase of 6-I3 is changed at 0.5 GPa to the massless 
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FIG. 4: (Color online) Angular dependence of the averaged 13 C-NMR shift d (central line shift of 
a Pake doublet) at the two unequal molecules in the unit cell (symbols) at (a) H/ jac (100 K) and 
(b) H/ /be (100 K). The same symbols are used as in Figs. 2 and 3. Curves stand for the least 
square fits to the data calculated from the crystal structure 3 - 3 - and the hyperfine-shift tensor given 
in the text. 

Dirac fermions with a linear energy-momentum dispersion, which is confirmed later by 13 C- 
NMR shift and relaxation-rate measurements performed at 0.8 GPa3. Noticeably, NMR 
spectra in the high-pressure phase exhibit more than eight lines. The qualitative difference 
in the spectral feature below and above 0.5 GPa evidences that the transition at 0.5 GPa 
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FIG. 5: (Color online) (a) Temperature dependence of the 13 C-NMR spectra under H//c in the 
2D conducting plane, (b), (c): Temperature dependence of the Knight shift, K (= 5 — a; with the 
central NMR line shift d and the temperature- independent chemical shift a ~ 106 ppm) and the 
nuclear spin-lattice relxation rate, 1/T\, divided by temperature T under the same field orientation 
(H//c). 



accompanies a structural phase transition with a symmetry reduction; thereby 8-l 3 is very 
probably transformed into a new structure similar to a-\ 3 that accommodates the massless 
Dirac fermions. 



B. Knight shift, nuclear spin-lattice relaxation rate, and electron-correlation effect 

Next, we turn our attentions to the local electronic structures in #-I 3 at ambient pressure. 
We measured the temperature dependence of 13 C-NMR spectra, which gives Knight shift, 
K, and nuclear spin-lattice relaxation rate, 1/T\. 

Figure |5](a) shows the temperature dependence of the 13 C-NMR spectra under an external 
field H of 6.00 T applied parallel to the c axis in the conduction layer. All molecules become 
equivalent under this field orientation, so only a single Pake doublet is observed. 

The spectra show moderate temperature dependence with decreasing temperature. The 
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NMR shift, S, defined as the midpoint of the doublet, is composed of the Knight shift, K, and 
the chemical shift, a, terms, i.e., 5 = K + a. The Knight shift, K, is related to the electron 
spin susceptibility Xs via K = aXs/2/iB^A, where a is the hyperfine-coupling constant, /^b is 
the Bohr magneton, and N\ is Avogadoro's number. The chemical shift, a, originates from 
the orbital motion of electrons in a BEDT-TTF molecule and is temperature independent. 
In the present field configuration, it is estimated as cr ~ 106 ppm from the chemical-shift 
tensors reported by Kawai et alM^- and the pseudo-orthorhombic crystal structure^ 3 -. By 
subtracting cr from the NMR shift 5, we evaluated the Knight shift K (= 5 — cr) as shown in 
Fig.E^b). It decreases monotonically with temperatures with a slight concave around 200 K, 
approaching a value of 100 ppm in the low-temperature limit. The temperature dependence 
of Knight shift K is more prominent than that of the spin susceptibility Xs determined by 
magnetization measurements earlier—. 

Figure[5]^c) shows the temperature dependence of nuclear spin-lattice relaxation rate 1/21 
divided by temperature, 1/21T, under the same field orientation (H//c). \jT{T exhibits 
a monotonic decrease with decreasing temperature. Compared with the Knight shift K 
shown in Fig. |5^b), 1/21 T shows relatively large temperature dependence, which seems to 
get steeper above ca. 200 K. In a metallic system, 1/Ti is proportional to the scattering 
rate of conducting electrons by nuclear spins near the Fermi level E-p and is known, in case 
of isotropic hyperfine interaction characterized by the coupling constant a, to follow the 
Korringa relation 38 — l/T x T = A7Tk B /h{-f n /-f c ) 2 a 2 {fi B /2N A ) 2 < D(E F ) 2 > T — where D(E F ) 
is the electronic density of states at Ep, 7 e is the gyromagnetic ratio of electron, and < x >t 
stands for the thermal average of the quantity x. In a conventional metal without electron- 
electron correlations, the following relation holds between Knight shift K and relaxation 
rate 1/21; (l/TiT)if- 2 (ft/47r£;B)(7e/7n) 2 = 1. In (BEDT-TTF) 2 X compounds, however, 
the hyperfine interaction is strongly anisotropic at the 13 C positions because of a significant 
contribution of dipolar interactions from the p z orbitals. This effect introduces a geometrical 
form factor in the above relation, namely, (l/TiT)K~ 2 (h/ '47r/cB)(7e/7n) 2 = (3{(, rj)^^ with 



^£ ^ (a xx /a zz ) 2 (sm 2 r] + cos 2 C cos 2 ^ + (a vy / a zz ) 2 (cos 2 r] + cos 2 C sin 2 //) + sin 2 ^ . 
2[(a xx /a zz ) sin 2 £ cos 2 rj + (a yy /a zz ) sin 2 £ sin 2 rj + cos 2 £] 2 



where (a xx , a yy , a zz ) are the principal values of the 13 C-hyperfine-coupling tensor, £ is the 
angle between the external field H and the z-principal axis shown in Fig. (He), and rj is 
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FIG. 6: (Color online) Temperature dependence of the 13 C-NMR Korringa ratio, K, oc (1/T\TK 2 ), 
under H/ jc. The dashed line stands for the free electron's case, K, = 1 (see the text). 

the polar angle measured from the x-principal axis in the xy plane of Fig. [U(c). (3{(,r)) is 
dependent on the orientation of the applied field H relative to the molecular axes. Further- 
more, the effect of the electron correlations is incorporated into the form by the introduction 
of the factor, IC, called the Korringa ratio which measures the degree of spin fluctuations 
quantitatively^, i.e., {l/TiT)K~ 2 {h/ '47rfcB)(7e/7n) 2 = /3(C) r i)^ ■ is unity for free electrons, 
while it becomes larger (smaller) than unity if there are anitiferromagnetic (ferromagnetic) 
spin fluctuations^. Using the hyperfine-shift tensor determined in Section III. A, the suscep- 
tibility Xs at 100 K reported by Salameh et al— (x s = 4.9 x 10~ 4 emu/mol formula unit), 
and the chemical-shift tensors determined by Kawai et al.—^-, we estimated the principal 
values of the 13 C-hyperfine-coupling tensor as (-1.1, -0.9, 5.6) (in kOe//iB)- Then, f3((,rj) is 
evaluated as ~ 0.97 with (C,^) = (139°, 85°) under the current field configuration (H//c). 

With the use of Knight shift K and \jT{T in Figs. Mjo) andE^c), and the /3((,r)) value 
determined above, the temperature dependence of the Korringa ratio /C is obtained as shown 
in Fig. [6j There is no strong temperature dependence in the Korringa ratio /C, which is 
2.4±0.2 up to 200 K and tends to increase into a value of 3.0 around room temperature. 
The value greater than unity points to the presence of antiferromagnetic spin fluctuations 
but is appreciably smaller than the typical values in strongly correlated electron systems 
realized in similar organic conductors. For instance, the 2D metal with strong electron 
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correlations, /t-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and K-(BEDT-TTF) 2 Cu(NCS) 2 i2 ^, locating 
in the vicinity of Mott transition 2 , show the Korringa ratio of K, ~ 8 at ambient pressure. 
The present value of /C hence indicates that electrons in #-I 3 are moderately correlated. This 
is consistent with the experiments of quantum oscillations^ - — giving an effective mass which 
is approximately a half of the K-type salts'. 

The overall results show no signature of charge ordering unlike most of the #-type com- 
pounds. However, it should be noted that a steep increase in \jT\T above 200 K [Fig. [5(c)] 
is a feature not seen in other metallic salts^and implies an enhancement of spin scatterings. 
The feature that the enhancement is less prominent in K, (Fig. [HD suggests that the spin 
scatterings are distributing over a wide range in the fe-space and are irrelevant to antiferro- 
magnetic instabilities. Recently, Cano-Cortes and coworkers^ theoretically investigated the 
extended Hubbard model with a modeled q-type structure. Among the consequences is the 
possible quantum phase transition of charge order with respect to the strength of electron 
correlations^ 1 ^. In the metallic side in the phase diagram, the system is predicted to exhibit 
a crossover from a Fermi liquid to an incoherent bad metal with quantum critical charge 
fluctuations at a finite temperature T*, which depends on the distance from the quantum 
critical point in the phase diagram. In the incoherent state above T*, the spin lifetime is 
expected to be shortened without preferential scattering vector. The turnabout behavior in 
1/TjT around 200K can be a signature of this crossover. Assuming that the largest transfer 
integral, t a , in 9-l 3 is 80 meV as deduced from the analysis of optical datar^, 200 K cor- 
responds to a reduced temperature of T*/t a = 0.2, which is in a reasonable range in the 
prediction 44 . A related feature is seen in the optical conductivity measurements by Take- 
naka et at— who observed a rapid loss of electronic coherence with increasing temperature. 
These experimental and theoretical results suggest that 9-1% is situated not so far away from 
the charge-ordering point. 



IV. CONCLUSION 

We performed 13 C-NMR measurements for a single crystal of #-(BEDT-TTF) 2 I 3 at am- 
bient pressure to probe microscopically the nature of the metallic state, which is a rare 
case in the #-type family of compounds. The orientation dependence of the NMR spec- 
tra revealed that the lattice symmetry can be well understood by the pseudo-orthorhombic 
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symmetry with a homogeneous molecular arrangement. The NMR shift and relaxation rate 
1/Ti measurements showed the Korringa relation up to room temperature with the weakly 
temperature-dependent Korringa ratio JC (= 2 ~ 3), which indicates that #-(BEDT-TTF) 2 i3 
is in a moderately correlated regime. However, the enhancement of relaxation rate, observed 
above 200 K, is a possible signature of the quantum critical charge fluctuations as suggested 
by a recent theoretical study. 
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